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We have explored the magnetic excitation spectrum of the S=l/2 square lattice Heisenberg an- 
tiferromagnet, K2V3O8 using both triple-axis and time-of-flight inelastic neutron scattering. The 
long-wavelength spin waves are consistent with the previously determined Hamiltonian for this ma- 
terial. A small energy gap of 72±9 fieV is observed at the antiferromagnetic zone center and the 
near-neighbor exchange constant is determined to be 1.08±0.03 meV. A finite ferromagnetic inter- 
planar coupling is observed along the crystallographic c-axis with a magnitude of J c =-0.0036±0.0006 
meV. However, upon approaching the zone boundary, the observed excitation spectrum deviates sig- 
nificantly from the expectation of linear spin wave theory resulting in split modes at the (ir/2,n/2) 
zone boundary point. The effects of magnon-phonon interaction, orbital degrees of freedom, multi- 
magnon scattering, and dilution/site randomness are considered in the context of the mode splitting. 
Unfortunately, no fully satisfactory explanation of this phenomenon is found and further theoretical 
and experimental work is needed. 



PACS numbers: 75.30.Ds, 75.50.Ee, 75.40.Gb 

I. INTRODUCTION 

Quantum magnetism has been a topic of considerable 
interest for many decades^, with particular interest in 
two-dimensional systems stimulated by the discovery of 
high-T c superconductivity in oxides comprised of Cu02 
layers 2 . Shortly after this discovery, it was realized that 
the parent compounds of these superconductors are well 
described by the quantum (S=l/2) square lattice Heisen- 
berg antiferromagnet (QSLHAF)^4 and since then this 
model system has been the topic of considerable theoret- 
ical and experimental investigation. 

Theoretically, there is growing consensus that the 
ground state of the QSLHAF is long-range ordered only 
at zero temperature. The dynamics of the QSLHAF 
are well described by classical linear spin-wave theory 
with quantum corrections, in the form of higher order 
expansions in necessary to extract accurate 



physical parameters. For instance, classical spin-wave 
theory accurately reproduces the dynamical structure 
factor, S(Q,w), with physical parameters, such as the 
spin-wave velocity, spin stiffness, and susceptibility 
strongly renormalized by these quantum corrections. 
Qualitative deviations from spin-wave theory may be 
manifest near the antiferromagnetic zone boundary, 
as shown from both quantum Monte-Carloi and series 
expansion-^ studies. These predictions suggest a 7-9% 
dispersion in the magnetic excitations between the 2d 
zone boundary points (tt/2,it/2) and (tt,0). Linear 
spin-wave theory predicts no dispersion between these 
points although recent extensions of the spin-wave 
calculation including up to second order corrections in 
1/S show a 2% dispersion^, much smaller than that 
predicted by other theoretical techniques. The effect 
of temperature on this model system is to destroy the 
T=0 long-range order and the temperature dependence 
of the correlation length, which provides a measure 
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of the degree of order, has been well studied both 
theoretically 10 ! 1 MVSand experimentally 14 ! 15 ! 16 ! 17 ! 18 . 
The temperature dependence of the dynamical structure 
factor has also been the topic of considerable theoretical 
effort 19 i 20 i 21 i 22 i 23 i 24 i 25 i 26 i 27 with much less experimental 
effort 28 ". In recent years theoretical investigations have 
been extended to include effects of an external applied 
magnetic field^ 9 *^ *^** 3 ^. Perhaps the most dramatic 
prediction of these studies is a decay of the single 
magnon spectrum into a two-magnon continuum for 
fields sufficiently close to the saturation field 31 . 

As is evident in the preceding paragraph, many of the 
interesting theoretical studies of the QSLHAF have in- 
volved detailed calculations of the dynamical properties. 
These properties are directly probed by inelastic neutron 
scattering experiments which measure the dynamical 
structure factor, S(Q,u/). Experimental measurements 
have focussed on the cuprates and inelastic neutron scat- 
tering measurements have been hindered by the large 
characteristic energy scale (J ~ 1500 K) of these mate- 
rials. Early reactor-based measurements on the high-T c 
parent compound, La2Cu04, were only able to measure 
the spin- wave velocity in the long-wavelength limils 4 * 3 * 3 * 3 -! 
The abundance of epithermal neutrons at spallation neu- 
tron sources have allowed for measurements of the full 
dispersion^ 5 ., resulting in a near-neighbor coupling con- 
stant of 136 meV for La2Cu04. Advancement in time-of- 
flight instrumentation and visualization software allowed 
much more detailed measurements to be performed 36 . 
In these experiments, the zone boundary dispersion was 
measured and a strong dispersion was observed between 
the (7r/2,7r/2) and (jrfi) zone boundary point but in the 
opposite trend to that predicted theoretically 6 * 7 * 8 . This 
trend was attributed to the presence of a higher order ring 
exchange interaction around CU4O4 square plaquettes 3 * 6 -, 
an interpretation which has been the topic of considerable 
interest both theoretically* 3 ^* 3 * 8 * 3 ** and experimentally 40 . 
In addition, several other related cuprate model systems 
have been studied. One of the best realizations of the 
QSLHAF is Sr2Cu02Cl2 where neutron scattering exper- 
iments have produced detailed measurements of the tem- 
perature dependent correlation length* 6 " but no studies 
of the dynamics have been performed to date. Another 
related compound, S^CuaC^Ck) is complicated by the 
presence of two interpenetrating square Cu 2+ sublattices 
with very different coupling constants (Jj — 130 meV 
and J 11 — 10 meV) and frustrated interactions between 
the sublatticesi 41 ". The magnetic excitation spectrum was 
measured in detail, for the specific range of temperatures 
where this material behaves as a QSLHAF with J ~ 10 
meV, and these measurements provided the first experi- 
mental evidence for the theoretically predicted dispersion 
between the (7r/2,7r/2) and (ir,0) zone boundary points 41 ". 

Temperature dependent measurements were limited by 
the large coupling constants in the cuprates to tempera- 
tures considerably less than J. These measurements were 
extended following the discovery of the metallo-organic 
Cu(DCOO) 2 -4D 2 (CFTD) which has a much smaller 



coupling constant, J — 6.3 meV— . Measurements of the 
correlation length were extended to temperatures com- 
parable to J 17 i 18 a nd were found to be in good agreement 
with theoretical predictions. Inelastic neutron scattering 
measurements on this material 2 *^ confirmed the theoret- 
ically predicted dispersion along the zone-boundary*?* 7 * 8 " 
and also measured the temperature dependence of the 
dynamical structure factor. Excitations were found to 
persist to temperatures as high as J/2 but were found to 
broaden and soften upon warming 28 . 

Experimental studies of the effects of an external mag- 
netic field have been impossible due to the magnitude of 
the coupling constants in known model materials. Inter- 
esting theoretical predictions in the presence of an ap- 
plied magnetic field require field strengths near satura- 
tion. Using the value of J, we can estimate the field re- 
quired for saturation for the model systems mentioned 
above. The resulting saturation field for La2CuC>4 is 
-3800 T, Sr 2 Cu30 4 Cl2 is -350 T (-4500 T) assuming 
the lower (higher) of the two coupling constants, and 
CFTD is —220 T. Clearly all of these materials are un- 
suitable for magnetic field measurements. Recently, a 
family of metallo-organics have been synthesized with 
much smaller coupling constants giving rise to saturation 
fields ranging from 2.25-25 Ti 2 .. Unfortunately, these ma- 
terials have rather strong three-dimensional interactions 
and are also problematic for inelastic neutron scatter- 
ing experiments as it has proven very difficult to grow 
large, deuterated single crystals. One particularly at- 
tractive material, and the topic of this communication, 
is K2V3O8 which has a saturation field of —38 T, exper- 
imentally achievable for some experiments. 

In this manuscript, we present detailed inelastic neu- 
tron scattering measurements and analysis of the zero- 
field magnetic excitation spectrum of K^VsOs- Such 
studies are not only essential to understand future re- 
sults obtained in the presence of a finite applied magnetic 
field but are quite interesting in their own right. Com- 
parison of the measured spectrum to the predictions of 
classical linear spin- wave theory has proven very interest- 
ing in other model systems as emphasized by the oppos- 
ing trend in the zone boundary dispersions observed in 
La2Cu04 36 and CFTD* 7 ". The zone boundary dispersion 
observed K2V3O8 is different than that observed in any 
of the previously explored model systems or predicted 
in theoretical studies, displaying two modes in the ex- 
citation spectrum in the immediate vicinity of the zone 
boundary. The remainder of the paper will be organized 
in the following manner; we begin by discussing the prop- 
erties of K2V3O8 (Section II) followed by a description of 
the details of the inelastic neutron scattering experiments 
(Section III). The experimental results are described in 
Section IV for spin-waves in the long-wavelength limit 
which is fully consistent with the prediction of linear spin- 
wave theory. In Section V the experimental results near 
the zone boundary are described showing the evidence for 
the two modes in the excitation spectrum. Several possi- 
ble sources of the observed zone boundary mode splitting 




FIG. 1: (Color online) (a) Crystal structure of K2V3O8 show- 
ing red V 4+ 05 pyramids linked by non-magnetic V 5+ 04 tetra- 
hedra shown in green. The grey spheres represent K + ions 
which separate the V3O8 2d sheets, (b) Projection of the 
crystal structure perpendicular to the c-axis. The intralayer 
coupling between S=l/2 V 4+ ions is shown by blue arrows, 
(c) Projection of the crystal structure perpendicular to the 
a-axis. The interlayer coupling, J c is shown by blue arrows. 



will be discussed in Section VI followed by conclusions in 
Section VII. 



II. K 2 V 3 8 : GENERAL PROPERTIES 

K2V3O8 crystallizes in the fresnoite structure with 
tetragonal space group P4bm and lattice constants 
a=8.87 A and c=5.2 A^. The room temperature struc- 
ture, shown in Fig. ^ consists of vanadium oxide planes 
composed of S=l/2, V 4+ -C>5 square pyramids and non- 
magnetic V 5+ -04 tetrahedra. The intralayer coupling 
constant, J, is small due to the long V 4+ -0-0-V 4+ su- 
perexchange path as can be seen in Fig. \T]p which also 
shows the square arrangement of the V 4+ ions. The 2d 
nature of the magnetic system results from layers of K + 
ions which separate the vanadium oxide layers. Pow- 
der magnetization measurements were well described by 
the QSLHAF with a small near neighbor coupling con- 
stant of 1.09 meV—. This material is found to order 
magnetically with Tjv=4 indicating a finite inter- 
layer coupling J c (Fig. ^). The magnetic structure 
in the absence of an applied magnetic field is a simple 
two-sublattice antiferromagnetic arrangement with spins 
aligned along the crystallographic c-axis with ferromag- 
netic stacking of the 2d layers^. A non-zero magnetic 
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FIG. 2: Reciprocal space diagram for K2V3O8. Structural 
Bragg reflections are indicated by black circles and mag- 
netic reflections by gray circles. Several high symmetry 
zone boundary points are indicated by gray diamonds. All 
points are labelled by the reciprocal space indices relevant to 
K2V3O8 (h and k). For reference, several points are addition- 
ally labelled using the square lattice notation. As can be seen 
from this notation, the square lattice is rotated by 45° from 
the lattice of K2V 3 Os. The dashed lines show the antiferro- 
magnetic zone boundary around the (1,0) zone center. 



field results in phase transitions when applied along both 
the c-axis and in the tetragonal basal plane. The o 
axis field induced phase transition is a rather conven- 
tional spin-flop transition while the presence of a basal 
plane applied field induces a peculiar spin reorientation 
transition where the spins rotate from the easy c-axis 
into the basal plane while remaining normal to the ap- 
plied field direction. These observations were described 
by the following two-dimensional Hamiltonian which in- 
cludes a Dzyaloshinskii-Moriya (DM) interaction and an 
additional c-axis anisotropy^, 

where represents near-neighbor pairs of spins, and 
S,, Sj are the spin operators for sites The applied 
field required to induce both the basal plane spin reori- 
entation and the c-axis spin flop transition allowed for 
estimates of E z /J of 0.0012 and D z /J of 0.04. 

Thermal conductivity measurements on K2V3O8 in- 
dicate dramatic enhancements in the low temperature 
heat transport in the presence of an applied magnetic 
field^. These same measurements indicate anomalies at 
110 K 4 4 suggesting a phase transition at this tempera- 
ture. Optical measurements clearly show anomalies in 
the local structure of the VO5 square pyramids on pass- 
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ing through 110 K 47 indicating the structural nature of 
this phase transition. 

Recently, a more comprehensive theoretical treatment 
of the magnetic properties of K2V3O8 has been carried 
outa£. This work pointed out that weak ferromagnetism 
induced by the DM interaction is "hidden" in K2V3O8 in 
the zero field limit by a dominant c-axis anisotropy. The 
effects of the DM interaction are revealed in the presence 
of an applied magnetic field and analytic expressions for 
the field strength required to induce both the spin-flop 
and spin reorientation transitions were provided. The 
field dependence of gap energies in the magnetic excita- 
tions were calculated and the effects of quantum correc- 
tions were included. One interesting theoretical predic- 
tion is a non-linear form for the gap energy as a function 
of field applied in the basal plane. 

The topic of this manuscript will be inelastic neutron 
scattering measurements in the limit of zero applied mag- 
netic field. If we set H=0 in Eq. ^and add a ferromag- 
netic coupling between the 2d layers, as expected from 
the known zero-field magnetic structure, linear spin- wave 
theory applied to the resulting Hamiltonian yields two 
degenerate spin wave modes with dispersion, 



wq = 2 J 
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where 



(2) 



7q,| = cos(fi,7r) cos(/c7r); 7q^ = sin 2 (/7r). (3) 



Here, we introduce an effective coupling constant, J = 
Z C J where Z c represents the quantum renormalization 
of the coupling constant with a best theoretical estimate 
of Z C =1.1S. 5 - 6 . Note that as a consistency check, if we set 
Jc=E z =D z =0, we recover the result for the QSLHAF, 



ujq = 2Jyl — cos 2 (/i7r) cos 2 (fc7r). 



(4) 



where the unit cell is rotated by 45° from the conven- 
tional square lattice unit cell (see Fig. |2J|. 

Inelastic neutron scattering experiments measure a 
convolution of the instrumental resolution function with 
the partial differential cross-section, 

d '° I F (Q) | 2 [n(u) + 1] (1+cos 2 <p)S(Q, u) (5) 



dfldoj 



where ki and fc/ are the magnitudes of the incident and 
final neutron wavevectors, F(Q) is the magnetic form fac- 
tor, in this case for V 4+ , n(u>) -I- 1 is the Bose occupation 
factor, (1 + cos 2 ip) is a polarization term reflecting the 
spin direction (<p is the angle between Q and the easy c- 
axis) and S(Q, ui) is the dynamical structure factor. Spin 
wave theory applied to the above Hamiltonian yields 
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(6) 




0.055 
0.05 
0.045 
0.04 
0.035 
[ 0.03 
0.025 
0.02 
0.015 
0.01 
0.005 



1.5 2 2.5 .5 1 

Q along (hOO) [r.l.u ] 



1,5 2 2.5 



FIG. 3: (Color online) Low energy magnetic excitation spec- 
trum of K2V3O8 at temperatures of 1.5K and 30K. The emer- 
gence of the excitations on cooling below the peak in the sus- 
ceptibility clearly shows the magnetic nature of the excita- 
tions. The temperature independent feature at H=1.8 and 
E=0.25 meV is spurious and is related to the presence of the 
nearby (2,0,1) structural Bragg reflection. The color bar in- 
dicates intensity of scattered neutrons in arbitrary units. 



III. EXPERIMENTAL DETAILS 

About 100 grams of V 2 5 (99.995%) and 76 grams of 
K2CO3 (99.997%) was loaded into a 250 ml Pt crucible 
and slowly heated in air to 700° C and held for 2 hours. 
The powder was added in two stages because of substan- 
tial foaming. The molten KVO3 flux was then cooled 
to room temperature and the Pt crucible and solidified 
KVO3 loaded into a fused silica ampoule that is necked 
down at the top. Using a long funnel, 9 grams of VO2 
(99%) was added to the Pt crucible and the entire am- 
poule was evacuated and sealed. The sealed ampoule was 
loaded into a furnace, heated to 850 °C for 6 hours, cooled 
to 700 C over lh and then cooled to about 400 °C at 1 
°C/h, followed by furnace cooling to room temperature. 
The KVO3 flux was then removed with a combination 
of warm water and ultrasonic vibration. The resulting 
K2V3O8 crystals are black rectangular plates with typi- 
cal dimensions of 5 x 5 x 1 mm 3 . 

Inelastic neutron scattering experiments on single crys- 
tal samples were performed using the backscattering 
spectrometer, IRIS at the ISIS neutron scattering facil- 
ity, the RITA II triple-axis spectrometer at the Paul- 
Scherrer Institut and the SPINS triple-axis spectrometer 
at the NIST Center for Neutron Research. The IRIS and 
RITA II experiments were performed on an array of ~ 
40 single crystals visually co-aligned with a total mosaic 
spread of about 4° . Improvements in crystal size allowed 
the SPINS experiment to be performed on an array of 
5 single crystals which were much more accurately co- 
aligned with a final mosaic spread of about 0.5°. For all 



5 




0.2 0.4 0.6 0.8 1.0 

Energy [meV] 

FIG. 4: (Color online) Cut along the energy axis for H ranging 
from 0.9 to 1.1 showing the clear presence of a small energy 
gap of about 70 fj,eV near the (1 1) antiferromagnetic zone 
center. The solid line represents the cut through simulated 
data resulting in extremely good quantitative agreement be- 
tween the data and the model. The strange lineshape is a 
result of the small number of detectors which cut through the 
dispersion at discrete positions. 



experiments, the crystals were mounted on an aluminum 
sample mount which was connected to the cold finger of 
a 4 He cryostat with a minimum temperature of about 
1.5 K. The IRIS experiments used PG (002) to select a 
fixed Ef of 1.845 meV yielding a resolution of ~18 /ieV 
at the elastic position. Measurements were performed 
using two sets of chopper settings allowing for studies 
with energy transfers covering a range from -0.2 to 1.2 
meV and 0.5 to 2.8 meV respectively. Both the RITA 
II and SPINS triple-axis spectrometers used PG (002) 
as the monochromator with pre-monochromator collima- 
tion determined by the presence of a neutron guide and 
an 80' post-monochromator collimator. The majority 
of measurements were performed with a fixed final en- 
ergy of 5 meV and a cooled Be filter before the analyzer 
resulting in an energy resolution of ^0.25 meV at the 
elastic position. Both instruments employ a multiblade 
PG (002) analyzer configuration and measurements were 
made with this analyzer array either in a flat or focussing 
configuration. 

For means of reference, we have included a map of 
reciprocal space in Figure with nuclear zone centers in- 
dicated by black circles, antiferromagnetic zone centers 
by gray circles, and several high symmetry antiferromag- 
netic zone boundary points by gray diamonds. To allow 
for easier comparison with previous measurements, the 
relevant region of reciprocal space is labelled by both h 
and k for K2V3O8 together with the labels for the con- 
ventional square lattice. This labelling clearly shows that 
the unit cell of K2V3O8 is rotated by 45° from the square 
lattice unit cell. 
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FIG. 5: (Color online) Measured excitation spectrum at sev- 
eral different sample rotations showing the dispersion along 
the crystallographic c-axis between the planes. The color 
bar represents the intensity of scattered neutrons in arbitrary 
units. 



IV. RESULTS: LONG WAVELENGTH 
EXCITATIONS 

One clear prediction of the dispersion shown in Eq. [21 
is the presence of an energy gap at the antiferromagnetic 
zone center of, 



A = 2 J 



\ 



1 + f' 



(7) 



Using the previously determined experimental values for 
J, E z /J, and D z /J^, together with the expected quan- 
tum correction J=1.18J£& we expect this energy gap 
to be about 73 /j,eV. In order to explore the low en- 
ergy region of the dispersion, we performed initial mea- 
surements on the IRIS backscattering spectrometer. For 
these measurements, the sample was mounted in the (h 
I) scattering plane. Measurements of the low energy 
portion of the magnetic excitation spectrum at tempera- 
tures of 1.5 K and 30 K are shown in Figure^] This plot 
shows the inelastic intensity as a color map as a function 
of component of momentum transfer along the (h 0) 
direction and energy transfer. The excitations, clearly 
seen in the T=1.5 K data, are found to vanish at higher 
temperatures providing clear evidence of their magnetic 
nature. The temperature independent bright intensity 
spot observed near h=1.8 and E=0.25 meV is a spurious 
feature, not uncommon in an inverse geometry, time-of- 
fllght spectrometer where a white beam is incident on 
the sample, related to a nearby (2,0,1) structural Bragg 
reflection. 

Careful examination of Figure [3] suggests a slight re- 
duction in inelastic intensity as the energy approaches 
the elastic position which may be an indication of a gap 
in the excitation spectrum. To see this more clearly, an 



6 




0.0 0.1 0.2 0.3 0.4 0.5 
Energy (meV) 



FIG. 6: (Color online) Comparison of the measured data at 
several different values of L with the simulated scattering. 
The gray bars represent the dispersion at the respective L 
values. This results in a ferromagnetic interplanar coupling 
of J c /J=-0.0028±0.0005 meV showing the nearly 2d nature 
of the magnetic system. 



energy dependent background was subtracted from the 
data shown in Figure |3| and an energy dependent cut 
through this background subtracted data for h ranging 
from 0.9 to 1.1 is shown in Figure This plot clearly 
shows a gap in the excitation spectrum which appears 
to about 100 fxeV. The solid line shown in Figure rep- 
resents the cross-section of Eq. [S] convoluted with the 
instrumental resolution function of the IRIS spectrome- 
ter. The line is in excellent agreement with the measured 
data and allows extraction of a zone center energy gap of 
72±9 fj,eV. As can be clearly seen from Eq. [7| the value 
of the gap cannot uniquely define both E z / J and D z j J. 
If we arbitrarily fix D z /J to value of 0.04^ and use a 
value of J of 1.275 (as will be shown later), we extract 
a value of the c-axis anisotropy of E z / J=0.0012±0.0001 
consistent with the value determined previously from the 
locations of the field-dependent phase transitions^. On 



a technical note, the apparent jumps in the data shown 
in Figure 0] result from the small number of detectors 
present on IRIS and these jumps are quite accurately re- 
produced in the simulated scattering, shown by the solid 
line in Figure 01 when the instrumental resolution and 
detector configuration are properly taken into account. 

The same experimental configuration on IRIS can also 
be used to measure the dispersion perpendicular to the 
2d planes. This is accomplished by measuring the same 
spectrum at several different sample rotations and the re- 
sults are shown for I w0.2,0.5 and 0.7 in Figure |SJ Note 
that these designations for I are only valid in the elas- 
tic position and I actually varies across the contour plot. 
Nonetheless, the dispersion along the c-axis can clearly be 
seen in this plot as a variation in the gap represented by 
the arrows in Figure^! To quantitatively extract the dis- 
persion along the c-axis, the same cut along energy (for h 
ranging from 0.9 to 1.1) at several different sample angles 
is shown in Figure EJ The solid lines shown in the fig- 
ure represent a numerical convolution of the cross-section 
(Eq. [SJ with the instrumental resolution. Once again, 
we see good agreement between the measured data and 
the simulation for all sample rotation angles measured. 
The grey vertical bars shown in each panel represent the 
location of the energy gap for the appropriate /-value 
clearly showing the c-axis dispersion. The fits allow us 
to extract a value of the interplanar coupling of J c / J=- 
0.0028±0.0005 reflecting the two dimensionality of the 
material. Note that the negative sign indicates the fer- 
romagnetic nature of the interaction consistent with the 
observed magnetic structure^. 

While IRIS is an ideal instrument for measuring fea- 
tures associated with the energy gap, due to extremely 
good energy resolution, it has Q-resolution limitations 
resulting from the small number of detectors on the in- 
strument and the large angle accepted by each analyzer 
crystal each of which scatters into a single detector. To 
see this, we plot cuts along the (hOO) direction for sev- 
eral ranges of energy transfer in Figure {7\ together with 
the corresponding numerical convolution of the expected 
cross-section (Eq. [5J| with the instrumental resolution 
function. While the simulation is seen to agree well with 
the measured data, it is extremely difficult to resolve the 
excitations and, hence, to extract physical parameters 
from these fits due to limited Q resolution. 

To extract information about the coupling constant, 
measurements were performed on the SPINS and RITA 
II triple-axis spectrometers. The multiblade analyzer 
system was configured in a horizontally focussing con- 
figuration for these measurements which increases inten- 
sity while maintaining good energy resolution at the sac- 
rifice of Q resolution in the direction perpendicular to 
k f . The results shown in Figure |S| represent constant-E 
scans performed with the sample mounted in the (/i,fc,0) 
plane using the SPINS triple-axis spectrometer. These 
measurements were performed along the (/i,l,0) direction 
which corresponds to the direction from (7t,7t) (for h=0) 
toward the antiferromagnetic zone boundary, (tt/2,tt/2) 
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FIG. 7: (Color online) Comparison of the data and simulated 
scattering at several ranges of energies. As before, the model 
can be seen to agree with the data up to 1 meV. The limited Q 
resolution of IRIS makes it difficult to resolve the dispersing 
modes. 
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FIG. 8: (Color online) Constant-E scans measured using the 
SPINS triple-axis spectrometer for energies ranging from 0.5 
meV to 1.5 meV together with fits to S(Q,w) convolved with 
the instrumental resolution function. 



(for /i=0.5). The much better Q resolution of the triple- 
axis spectrometer can be seen from the clearly disper- 
sive mode in the data shown in Fig. |H1 The use of the 
horizontally focussing configuration made analysis of the 
data complicated in this case. Typically, this configura- 
tion is used with an irrelevant direction in the scattering 
plane which can be oriented along the direction perpen- 
dicular to kf. However, in these measurements, both 
directions within the (h,k,0) plane are strongly disper- 
sive. This complicated the analysis as typical Gaussian 
approximations to the resolution function were not able 
to reproduce the sharp features observed in the data. 
Consequently, the cross-section appropriate to K2V3O8 
(Eq. EJ), was convolved with the resolution function for 
each of the 7 individual analyzer blades and the contribu- 
tion from each blade was added with the correct phase to 
yield the solid lines shown in Figure |H1 To correctly du- 
plicate the effective collimation introduced by the small 
width of the individual analyzer blades, the more rigor- 
ous Popovici approximation^ to the triple-axis resolu- 
tion function was required. Several scans ranging from 
0.5 meV to 1.5 meV were simultaneously fit with a single 
overall amplitude clearly showing that intensity modula- 



tion represented by <S(Q, to) (Eq. 0) provides an excellent 
description of the data. This analysis allowed extraction 
of the 2d coupling constant J of 1.275±0.03 meV. 

As mentioned previously, the coupling constant, J, is 
written as an effective coupling and theoretical studies 
indicate that quantum corrections to spin-wave theory 
result in an overall renormalization of the dispersion by 
a scaling factor Z c such that J=Z c j£&. As a consistency 
check, we can calculate the value of Z c by comparing the 
coupling constant above to the value of J=1.09 meV ob- 
tained from magnetization measurements^* where quan- 
tum effects were taken into account. This results in 
Z c =1.174±0.03, in excellent agreement with the best the- 
oretical estimate of Z c of 1.18££ justifying the designa- 
tion of K2V3O8 as a good example of a QSLHAF. Using 
the theoretical value for Z c , we can extract the true cou- 
pling constant, J, for K 2 V 3 8 of 1.08±0.03 meV. 

Measurements of the dispersion were additionally per- 
formed along the square lattice direction (7r,7r) to (7r,0) 
corresponding to the (h,l+h,0) direction in K^VsOs- 
Measurements along this direction were performed us- 
ing the RITA II triple-axis spectrometer with the sample 
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FIG. 9: (Color online) Summary of measurements performed 
along the (h,l+h,0) direction using the RITA II triple-axis 
spectrometer is shown in the upper panel. The lower panel 
shows the result of the best fit simulation to S(Q,ai) convolved 
with the instrumental resolution function. The color bar rep- 
resents the intensity of scattered neutrons in counts per ~2 
minutes. 



mounted in the (h,k,0) scattering plane and the results 
are shown in the upper panel of Figure The solid line 
in the upper panel represents the classical spin wave dis- 
persion with J=1.275 meV which clearly describes the 
data well. The lower panel represents the best fit to a 
convolution of Eq. [S] with the instrumental resolution 
function with the dispersion represented by the solid line 
in the upper panel. The fits agree very well with the data 
across the entire zone except in the immediate vicinity 
of the zone boundary where the peak intensity appears 
weaker than the calculated peak value. The behavior 
near the zone boundary will be discussed in more detail 
in the next section. 



V. RESULTS: SHORT WAVELENGTH 
BEHAVIOR 

As mentioned in the introduction, the dynamical struc- 
ture factor of the QSLHAF only deviates from predic- 
tions of spin- wave theory near the antiferromagnetic zone 
boundary. Consequently, to explore the behavior in the 
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FIG. 10: (Color online) The measured dispersion from energy 
transfers ranging from 0.5 to 2.8 meV using the IRIS backscat- 
tering spectrometer is shown in the upper panel. The lower 
panel represents the results of a convolution of the instru- 
mental resolution function with the expected cross-section. 
In both panels, the identical classical spin-wave dispersion 
dispersion, corresponding to a J=1.275 meV is superposed. 
The color bar represents the intensity of scattered neutrons 
in arbitrary units. 



vicinity of the zone boundary more carefully, measure- 
ments on IRIS were extended with the choppers config- 
ured to select a range of energy transfers from 0.5 to 
2.8 meV and the resulting data are shown in the up- 
per panel of Figure As can be clearly seen from the 
data, we observe well defined excitations up to the zone 
boundary. However, as can be seen from the solid line 
in the upper panel, the dispersion curve which describes 
the data well near the zone center seems to disagree with 
the data near the zone boundary. In the lower panel of 
Figure we show the result of a convolution of the 
instrumental resolution with the expected cross-section 
together with the dispersion curve used to generate the 
simulation. In comparing this to the measured data, it 
is also clear that once the amplitude is set to provide 
good agreement with the data near the zone center, the 
calculated intensity near the zone boundary exceeds that 
of the measurement. However, closer examination of the 
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FIG. 12: (Color online) Summary of the full triple-axis mea- 
sured dispersion. The measurements are the result of experi- 
ments performed at both RITA II and SPINS for experiments 
in both the (h,k,0) and (h,0,l) scattering planes. The solid 
line corresponds to a fit to linear spin-wave theory for data 
near the zone center. The dashed line represents the quantum 
corrections to the dispersion—. 



FIG. 11: (Color online) Cut through the IRIS data for a range 
of h from 0.4 to 0.6 as a function of energy transfer with a 
nominal Q of (0.5,0,1.1) is shown in the upper panel. The 
lower panel shows the corresponding constant-Q scan through 
(0.5,0,1) using the RITA II triple-axis. The solid gray bar 
shows the range of zone boundary energies resulting from 
measurements closer to the zone center and the predictions 
of linear spin- wave theory. 



measured intensity indicates the presence of two modes 
with intensity in the second mode only appearing near 
the zone boundary. 

Instrumcntally, the IRIS spectrometer exhibits various 
spurious scattering features around 3 meV energy trans- 
fer. The presence of these features led to scepticism as 
to the validity of the observed two modes near the zone 
boundary. To confirm this observation, measurements 
were performed in the identical (hfl,l) scattering plane 
using the RITA II triple-axis spectrometer. The instru- 
ment was configured with E/=4.6 meV and the resulting 
constant-Q scan at Q=(0.5,0,l) is shown in the lower 
panel of Figure 1111 In the upper panel, we show a cut 
through the IRIS data from h=0A to h=0.6 with the sam- 
ple slightly rotated from the measurement shown in Fig- 
urern]lso as to emphasize the data near h=0.5. Both the 
time-of-flight and triple-axis data clearly show the pres- 
ence of two peaks near the zone boundary. For reference, 
the gray box represents the predicted zone boundary en- 
ergy corresponding to J= 1.275 meV with the width of 
the box corresponding to the uncertainty of ±0.03 meV. 
Clearly, in both measurements, the calculated value of 
the zone boundary energy falls between the two mea- 



sured peaks. Note that the exact energies of the two 
modes differ slightly between the two instruments, a dis- 
crepancy which can be accounted for by differences in 
energy calibration. 

The full dispersion, measured using triple-axis spec- 
trometers along high symmetry directions in the 2d 
plane, is summarized in Figure IT21 with Q plotted m 
square lattice notation. The data shown in this figure are 
the result of experiments performed using both the RITA 
II and SPINS triple-axis spectrometers with the sample 
oriented in both the (h,k,0) and (h,0,l) scattering planes. 
The measured data along the (tt,tt) - (7r,0) direction and 
near the zone center along the (tt,tt) - (ir/2,ir/2) were fit 
to the expected S(Q,cl>) convolved with the instrumental 
resolution function. The amplitude, coupling constant, 
and a Lorentzian broadening term were allowed to vary 
for each of these data points. As this model no longer 
describes the data well near the (tt/2,tt/2) zone bound- 
ary point, the data here was fit to Gaussian lineshapes 
with widths constrained to be broader than that expected 
from resolution. The solid line is the expected classi- 
cal dispersion (Eq. |2$ with J=1.275 meV, D 2 /J=0.04, 
E/J=0.0012, and J c / J=-0.0028. This solid line is seen 
to agree well with the data near the zone center but devi- 
ates on approaching the (n/2,ir/2) zone boundary point. 
As the data begins to deviate from the solid line, an ad- 
ditional mode appears. As reflected by the size of the 
error bars, this mode is rather weak when it first ap- 
pears and gains intensity on approaching the zone bound- 
ary. As shown previously in Figure 1111 the calculated 
dispersion falls between the two modes on approaching 
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(7r/2,7r/2). There is no clear evidence for two peaks near 
the (7T,0) zone boundary point. However, as indicated 
by the large error bars near this point, the measured 
spectrum is broader than would be expected from reso- 
lution. The dashed line represents the expected quantum 
dispersion^. While it does appear that the data agrees 
better with the classical dispersion, there very well may 
be two modes near the (tt,0) zone boundary point as well. 
As such, the points shown on the dispersion may repre- 
sent the mid point of the two modes the lower of which 
would be quite close to the quantum dispersion. 



As mentioned in the introduction, the one region where 
measurements and quantum calculations on the QSL- 
HAF deviate from the predictions of spin-wave theory 
is in measuring the zone boundary energy along the di- 
rection (7r/2,7r/2) to (7r,0). Consequently, we have mea- 
sured the excitations along this direction using the SPINS 
triple-axis spectrometer to study the evolution of the 
observed two modes and the result is shown in Figure 
IT31 For reference, as shown in Figure [21 Q=(0.5,l,0) 
corresponds to the square lattice point (ir/2,ir/2) while 
Q=(0.5,1.5,0) corresponds to (nfl). The presence of two 
modes can also be seen in this scattering plane. To fur- 
ther emphasize this, the solid lines in the figure are the 
results of two Gaussians centered at the positions mea- 
sured under higher resolution conditions, shown in Fig. 
1111 with fixed widths and only the amplitudes of the 
Gaussians varying. The results clearly show an inten- 
sity distribution which varies across the zone boundary 
with the higher energy mode having its greatest intensity 
near (0.5,1,0) and continuously decreasing on approach- 
ing (0.5,1.5,0). In fact, although this two mode model 
still describes the data fairly well near (0.5,1.5,0), a bet- 
ter description, as shown by the solid gray line is that of 
a broadened Gaussian. This is consistent with the pre- 
viously presented data showing the dispersion along the 
(h,l+h,0) direction (Fig. l§land lT2|) which suggests only 
a single mode along this direction albeit with increased 
width near the zone boundary. As mentioned previously, 
there may be two modes in this direction as well and 
the weaker intensity of the upper mode may act to in- 
crease the difficultly in resolving the individual modes. 
For reference, the solid gray bar in Figure IT31 shows the 
range of values for the zone boundary energy given by 
the dispersion of Eq. [SJwith J=1.275 meV. The width of 
this bar is reflective of the uncertainty in the predicted 
classical zone boundary energy and is approximately 5%. 
Theoretically, dispersion between (7r/2,7r/2) and (7r,0) is 
predicted to be 7-9%, less than twice the width of the 
bar. Clearly, the measured data near the zone boundary 
is far too broad to allow observation of such a small zone 
boundary dispersion and the presence of the unexpected 
two modes makes such an investigation meaningless. 
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FIG. 13: (Color online) Constant-Q scans along the zone 
boundary for Q values ranging from (0.5,1,0) to (0.5,1.5,0). 
The solid lines present fits to two Gaussians with positions 
fixed at values measured in higher resolution conditions with 
fixed widths and only the amplitudes varying at different Q 
positions. The gray solid line through the black data shows a 
fit to a single, broadened Gaussian. The solid gray bar shows 
the range of values of zone boundary energies predicted using 
spin-wave theory and consistent with measurements near the 
zone center. 



VI. ORIGIN OF SPLITTING 

A. Magnon-Phonon Interaction 

The observation of two modes near the zone bound- 
ary is clearly not expected from the Hamiltonian pos- 
tulated for K2V3O8 (Eq. P). We have experimentally 
characterized the zone boundary properties quite exten- 
sively. The higher energy mode only seems to have sig- 
nificant spectral weight in the immediate vicinity of the 
zone boundary. We have tried fitting the data closer to 
the zone center using a model incorporating two disper- 
sive modes and no significant improvement in fit quality 
was observed supporting the notion that the high en- 
ergy mode is only significant near the zone boundary. In 
addition, the presence of two modes is much clearer near 
Q=(0.5,l,0) (or (71-/2, 7r/2) in square lattice notation) and 
the intensity of the upper mode falls off continuously in 
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moving away from this position along the zone boundary. 
Eventually, this evolves into what is best described as a 
broadened single peak at Q=(0.5,1.5,0) ((7r,0) in square 
lattice notation). In addition, both modes vanish upon 
increasing the temperature to ^30 K indicating that both 
are of magnetic origin. This rules out the possibility of a 
low-lying phonon mode interacting with the spin waves. 
Any proposed model to describe the observed scattering 
must satisfy these experimental observations. 



B. Additional Magnetic Interactions 

Although we haven't determined a single, consistent 
explanation for the peculiar zone boundary properties 
of K2V3O8, there are a number of possible explanations 
which we can discuss. Of course, the first possibility is 
that of additional couplings in the Hamiltonian. Longer 
range interaction do, in fact, tend to manifest them- 
selves near the zone boundary 36 . Although we can't 
rule out such a possibility, it is difficult to conceive of 
a coupling term which will only produce intensity near 
the zone boundary. The particular case of next-near- 
neighbor coupling has been studied quite extensively us- 
ing both linear spin-wave theory^ and series expansion 
techniques^. Although the inclusion of such a coupling 
term does modify the dispersion, with the largest effects 
observed near the zone boundary no evidence for a split 
excitation spectrum is predicted£2*£i. The symmetry ar- 
guments presented by Bogdanov et alSL suggest that 
K2V3O8 together with Ba 2 CuGe20 , ^ 3 '^ 4 are members of 
an interesting class of compounds where chiral inhomo- 
geneous magnetic structures can coexist with weak fcr- 
romagnetism. The net result of this symmetry analysis 
is the existence of an additional, symmetry allowed term 
in the Hamiltonian describing K^VsOs- The form of this 
chiral coupling term is rather complicated and its effects 
on the excitation spectrum are unclear at this point in- 
dicating the need for further theoretical investigations. 



C. Multimagnon Scattering 

Another possible explanation for the additional scat- 
tering observed at the zone boundary is multi-magnon 
scattering. A recent theoretical investigation 3 — based on 
the assumption that high-energy fermions exist in the 
QSLHAF suggest a very structured multi-magnon con- 
tinuum resulting in a rather sharply peaked extra mode 
around the (7r/2,7r/2) point and much broader scattering 
around (7r,0) in qualitative agreement with our experi- 
mental observations. However, the strength of the scat- 
tering in the additional mode makes such an interpreta- 
tion very unlikely. Theoretical studies using both series 
expansion^ and quantum Monte Carloi indicate a sig- 
nificant amount of spectral weight in the multi-magnon 
spectrum. The calculated ratio of the longitudinal to 
transverse structure factors for unpolarized neutrons is 



29% (34%) at (tt,0) while at (W/2,7r/2) the ratio is 31% 
(35%) where the values quoted are the results of series 
expansion^ (quantum Monte Carloi) studies. Although 
these values are rather large reflecting the strength of the 
multi-magnon scattering expected for the QSLHAF, they 
are somewhat deceiving in that they assume that neutron 
scattering will fully measure two transverse modes and 
one longitudinal mode. In reality, the neutron only sees 
components of magnetic moment normal to the wavevec- 
tor transfer resulting in a correction term which has the 
form l+cos 2 (<p) for the two transverse modes and sin 2 (<£>) 
for the longitudinal mode where <p is the angle between 
Q and the spin direction, in this case the c-axis. For 
the results presented in Figure ITT1 the angle between the 
(0.5,0,1) direction and the c-axis is about 30 degrees re- 
sulting in a correction term of 1.75 for the two transverse 
modes and 0.25 for the longitudinal mode. This implies, 
at this wavevector, that instead of seeing about 30% of 
the intensity, as predicted for the calculation under the 
assumption of a 1:2 ratio of longitudinal to transverse 
we would actually see 30%*2*0.25/1.75 « 8.5%. Clearly, 
the measurement at (0.5,0,1) shows intensities for the 
two modes which are very similar in magnitude and is 
not consistent with the upper mode being 8.5% of the 
lower mode. This essentially eliminates multi-magnon 
scattering as a source for the additional mode. 



D. Orbital Effects 

Another possible source for the observed mode split- 
ting is the influence of orbital degrees of freedom. Ex- 
amination of the dispersion curve (Fig. 1121 is qualita- 
tively consistent with the behavior expected if a flat mode 
were interacting with the spin wave dispersion resulting 
in repulsion of the modes. In addition, a flat excitation 
with a weak neutron structure factor which was enhanced 
through mixing with the spin-wave mode could explain 
the observed intensity modulation. One potential source 
of such a flat mode would be a low-lying crystal field exci- 
tation. A related effect was predicted theoretically using 
a generalized Holstein Primakoff transformation of the 
magnetic excitations in a system with spin and orbital de- 
grees of freedom in the presence of orbital degeneracy^.. 
The net result of inclusion of orbital degrees of freedom 
is an excitation spectrum with modes of pure spin, pure 
orbital, and mixed spin-orbital character—. Mixing is 
predicted between the spin and spin-orbital modes the 
strength of which is greatest near the zone boundary 
and weaker near the zone center. The scenarios postu- 
lated above require an understanding of the crystal-field 
ground state together with the approximate energies for 
the lowest excitations from this ground state. 

Symmetry analysis of the local square-pyramidal envi- 
ronment of the V 4+ ion yields a splitting of the octahedral 
t2 g levels into a singlet (d xy ) and a doublet (d xz , d yz ) and 
the octahedral e g states into two singlets (d x 2_ y 2 and 
d z i). Consequently, the only possibility of an orbitally 
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degenerate ground state would the d xz , d yz doublet. In 
an ideal square pyramidal geometry, where the cation 
sits in the center of the square plane of the pyramid, 
this doublet is indeed the ground state. However, for 
V A+ 0§ square pyramids, the V A+ cation is significantly 
displaced from this square plane which lowers the energy 
of the d xy singlet relative to that of the doublet. Conse- 
quently, the detailed structural arrangement needs to be 
taken into account to determine the orbital ground state. 
Ohama et. a& calculated the crystal field splitting as 
a function of displacement of the cation from the basal 
plane under the assumption of a ideal square pyramid 
with equal V-0 distances of 1.9 A. This calculation sug- 
gests that when the distance exceeds about 0.35 A, the 
d xy singlet is the crystal field ground state. As the true 
distance for K2V3O8 is 0.55 A, this calculation suggests 
a non-degenerate singlet crystal field ground state ruling 
out the above scenario of a mixed spin-orbital mode. 

To rule out the possibility of low-lying crystal field 
modes, we must consider the splitting between the d xy 
ground state and the d xz , d yz lowest excited states in a 
more quantitative sense. Recent interest in S=l/2 sys- 
tems where quantum effects are significant have resulted 
in several materials exhibiting V 4+ 05 square pyramidal 
local structure^,. Detailed molecular orbital calculations 
for a'-NaV 2 5 , CaV 2 5 , and MgV 2 5 58 all of which ex- 
hibit square pyramidal geometry similar to K2V3O8 con- 
firms the d xy ground state and estimates a rather large 
splitting to the lowest excited state ranging from 600-900 
mcV. Evidence of this large crystal field splitting is also 
seen in band structure calculations for several square- 
pyramidal V 4+ vanadates where the d xy bands are well 
separated from the remainder of the d bands£2i£2i§ii§2iS. 
Similar calculations performed on K^VsOg^ show a very 
similar splitting of the d xy bands suggesting a rather large 
crystal field splitting. One particular compound of note 
is Li 2 VOSi04§£i which has a crystal structure where V-0 
distances within the square pyramid are almost identi- 
cal to K 2 V30s- For this compound, an experimental 
estimate of the splitting between the ground state and 
first crystal field excited state yielded a splitting of ~150 
meV 65 while LDA band structure calculations estimated 
this value to be ^400 meV— . The local structural sim- 
ilarity between K 2 V30s and these other vanadates to- 
gether with the evidence for rather large splitting in all 
of these materials, is very strong evidence of splitting 
between the d xy ground state and the d xz ,d yz lowest ly- 
ing excited state for K 2 V30s which is on the order of 
hundreds of meV. A splitting of this magnitude together 
with the non-degenerate orbital ground state eliminates 
orbital effects from playing a significant role on the mode 
splitting observed at much lower energies (~3 meV). As 
a caveat, it is important to note that the crystal field 
splitting in K 2 V30s has not been directly measured ex- 
perimentally and the above arguments are based purely 
on comparison with other related vanadate systems to- 
gether with band structure calculations. 



E. Dilution/Randomness 

The measured anomalies near the zone boundary in the 
magnetic excitation spectrum of K 2 V3 0s are strikingly 
similar to the calculated excitation spectrum for site di- 
luted 2d antiferromagnets^ 7 ^. Studies of the effects of 
disorder on lattice models^ showed the development of 
localized modes, named fractonsfii with increasing disor- 
der. Similar localized modes were predicted to exist in 
studies of the effects of disorder on Heisenberg ferromag- 
nets and antiferromagnets 69 . Fractons have little impact 
on the excitation spectrum at longer wavelengths as this 
effectively averages over enough sites to make dilution ef- 
fects negligible. However, as the characteristic distance 
becomes smaller on approaching the zone boundary, the 
influence of these fractons becomes much more signifi- 
cant and theoretical predictions 7 ^ suggest a multi-peaked 
spectrum on approaching the zone boundary. Experi- 
mentally, inelastic neutron scattering experiments on the 
3d antiferromagnet, Mn 2; Zni_ 2 ;F 2 S and the 2d antifer- 
romagnet, Rb 2 Mn T Mgi_ x F4H, provided clear evidence 
for the multi-peaked nature of S(Q,w) in the presence 
of substantial disorder with concentrations in both cases 
very close to the percolation limit. 

Recently, detailed calculations of S(Q,w) have 
been presented for the site diluted 2d Heisenberg 
antiferromagnetSLS. These calculations show behav- 
ior which is similar to the undiluted case at long wave- 
lengths, albeit with a renormalized spin- wave velocity^ 7 -, 
with a crossover on approaching the zone boundary to 
a multi-peaked spectrumSLS the form of which is very 
similar to the experimental observations presented here 
for K 2 V30g. Specifically, at the zone boundary, the the- 
oretical predictions show two modes at weaker doping 
concentrations, with a splitting between the modes which 
is largely independent of doping level. This theoretical 
splitting is ^0.6 meV using the experimental parame- 
ters of K 2 V3 08 which can be compared to the exper- 
imentally measured splitting in K 2 V30s of 0.45 meV. 
These values are in fairly good agreement, particularly 
given that anisotropy, interplanar coupling, and quan- 
tum effects haven't been accounted for in the theoretical 
models 67 ' 68 . There are, however, some differences be- 
tween the experimental observations and the theoretical 
predictions which should be pointed out. First, the cal- 
culated S(Q,lj)£L§£ contains two modes on approaching 
both the (7r/2, 7r/2) and (tv,Q) zone boundary points while 
the experimental observations only show a clear splitting 
near (ir/2,w/2) with what is best described as broadening 
near (7r,0). In addition, the calculations shown to date 
make no predictions of how S(Q,oj) varies along the zone 
boundary direction (7r/2,7r/2)-(7r,0). Consequently, it is 
unclear as to whether the experimentally observed spec- 
tral weight shift can be described by this model. Clearly, 
further theoretical work is needed to determine if the 
qualitative comparison to the site dilution model is able 
to describe the anomalous features seen in the excitation 
spectrum of K 2 V30s- 
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The obvious question that arises, however, is why there 
would be significant disorder in K2V3O8. This material 
is a mixed valence material with an average valence of 
4.667, containing both magnetic V 4+ and non-magnetic 
V 5+ . Crystallographically, the three Vanadium positions 
in each formula unit are comprised of 2 tetrahedrally co- 
ordinated sites and 1 pyramidally coordinated site. The 
simplest manner of getting the correct overall valence 
is to have the V 5+ located on the tetrahedral site and 
the V 4+ on the pyramidal site and this is the assumed 
charge ordering scheme. However, there is no strict ex- 
perimental proof that such an ordering occurs in K2V3O8 
and some valence mixing on the pyramidal site could 
result in disorder effects. Other mixed valence vana- 
dates, most notably a'-NaV205 are strongly influenced 
by charge disproportionation and the resulting charge de- 
grees of freedomS 1 * 7 ^. More detailed experimental studies 
are necessary to explore the possibility of site disorder in 
the low temperature state of K 2 V 3 8 . 



F. Structural Distortion 

One final experimental note of potential importance 
is the recent observation of a higher temperature phase 
transition (T c -110 K) in K 2 V 3 0fA This tran- 
sition is clearly structural in nature, as evidenced 
by changes in the vibrational properties— Com- 
pounds such as K2V3O8 with the fresnoite structure 
are typically sensitive to displacive structural phase 
transitions resulting in incommensurately modulated 
structures— In most systems, the incommensurate 
modulation is along the (1,1,0) crystallographic direc- 
tion and usually results in a modulation wavevector 
of Q~0.3(l,l,0) or Q~0.3(l,l,0)+0.5c*. Very recent 
crystallographic measurements on K2V3O8 indicates the 
presence of such a displacive transition with superlat- 
tice peaks seen in x-ray diffraction studies with wavevec- 
tor Q— 0.3(1, l,0)+0. 5c — . Experimentally, no direct ev- 
idence for this changed periodicity is observed in the 
inelastic neutron scattering measurements. The pres- 
ence of an incommensurate structural modulation will 
affect the interactions between near-neighbor spins via a 
long-wavelength modulation of the superexchange inter- 
actions. The effects of periodic modulations of exchange 
interactions has received some attention in recent years 
in investigating the effects of stripe order on the magnetic 
interactions in cuprates2S*2L2&. A recent theoretical in- 
vestigation described the modification of the Hamiltonian 
required to account for such an exchange constants mod- 
ulation in the limit of weak modulation strength 79 . This 
study suggests that simple model systems, for instance, 
conventional antiferromagnets such as K2V3O8, would be 
insensitive to small long-wavelength modulations of the 
exchange constants. It should be noted, however, that 
this study focussed on the ground state properties and 
not more subtle changes to the spin-wave spectrum. In 
addition, this study assumes a weak modulation of the 



exchange which need not be the case. More theoretical 
effort is needed to estimate the strength of the expected 
modulation for the particular case of K2V3O8 and the 
resulting effect of such a modulation on the spin-wave 
spectrum. 

Interestingly, the material RD2V3O8 which is isostruc- 
tural at room temperature, also exhibits a structural 
phase transition with a very different modulation of 
Q— 0.16c*— . These distortions are explained, within the 
rigid unit mode analysis, as resulting from rotations of 
the V 5+ 04 tetrahedra such as to uniformly raise or lower 
the V 4+ 05 pyramidal network along the c-axis22. The 
proposed modulation does not affect the V 4+ -0-0-V 4+ 
superexchange pathway within the 2d planes only affect- 
ing the separation between layers and, hence, the inter- 
planar coupling. One would expect such a distortion 
to have little effect on the excitation spectrum which is 
predominately two-dimensional in character. Compari- 
son of the magnetic excitations in Rf^VjOs to those of 
K2V3O8, where the modulation contains a large in-plane 
component, could shed some light on the nature of the 
observed mode splitting. 



VII. CONCLUSIONS 

In conclusion, we have performed detailed inelas- 
tic neutron scattering studies of the magnetic excita- 
tion spectrum of the 2d antiferromagnet, K2V3O8 us- 
ing a combination of triple-axis and time-of-flight exper- 
iments. The long-wavelength region of the excitation 
spectrum confirms that K2V3O8 is an excellent exam- 
ple of a QSLHAF with a very small near-neighbor cou- 
pling constant J=1.08±0.03 meV. In addition, we were 
able to confirm the presence of a very small anisotropy 
gap in the excitation spectrum of 72±9 fieV. Under 
the assumption of the previously suggested Hamilto- 
nian with the Dzyaloshinskii-Moriya interaction set to 
D z /J=0.04, we determine a c-axis exchange anisotropy 
of E/, 7=0.001 2 ±0.0001 in very good agreement with 
the values estimated based on the location of the field- 
induced phase transitions 4 ^.. Finally, the dispersion along 
the c-axis was measured indicating a ferromagnetic inter- 
planar coupling J c =-0.0036±0.0006 meV (using J c /J=- 
0.0028±0.0005 and J=1.275±0.03) demonstrating that 
K2V3O8 is a very good two-dimensional material. 

As we approach the zone boundary, however, the pre- 
viously determined model no longer adequately describes 
the measured excitation spectrum. The excitations near 
the (7r,0) zone boundary point are broadened while near 
the (7r/2,7r/2) zone boundary point we observe two clear 
modes rather evenly split around the expected posi- 
tion of the single mode based on measurements in long- 
wavelength limit and the predictions of spin-wave the- 
ory. The upper of the two modes has a peculiar Q- 
dependence in that it only has significant intensity near 
the zone boundary. In addition, the intensity of the up- 



14 



per mode is seen to fall off continuously as we move along 
the zone boundary from (ir/2,n/2) to (ir,0). We have dis- 
cussed several possible explanations for the split mode 
near the zone boundary. Among these explanations, the 
best agreement with the data is the case of disorder on 
the 2d lattice although we have no other evidence for 
disorder in this material. We hope that these measure- 
ments will stimulate additional effort both experimen- 
tally and theoretically to attempt to explain the peculiar 
zone boundary properties of K2V3O8. 
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